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ABSTRACT 



By means of a series of N-body/SPH simulations we investigate the morphological properties of thick stellar disks formed through 
minor mergers with, e.g. a range of gas-to-stellar mass ratios. We show that the vertical surface density profile of the post-merger 
thick disk follows a sech function and has an excess in the regions furthest away from the disk mid-plane (z > 2 kpc). This stellar 
excess also follows a sech function with a larger scale height than the main thick disk component (except at large radii). It is usually 
dominated by stars from the primary galaxy, but this depends on the orbital configuration. Stars in the excess have a rotational velocity 
lower than that of stars in the thick disk, and they may thus be confused with stars in the inner galactic halo, which can have a similar 
lag. Confirming previous results, the thick disk scale height increases with radius and the rate of its increase is smaller for more gas 
rich primary galaxies. On the contrary, the scale height of the stellar excess is independent of both radius and gas fraction. We also 
find that the post-merger thick disk has a radial scale length which is 10 - 50% larger than that of the thin disk. Two consecutive 
mergers have basically the same effect on heating the stellar disk as a single merger of the same total mass, i.e., the disk heating effect 
of a few consecutive mergers does not saturate but is cumulative. To investigate how thick disks produced through secular processes 
may differ from those produced by minor mergers, we also simulated gravitationally unstable gas-rich disks ("clumpy disks"). These 
clumpy disks do not produce either a stellar excess or a ratio of thick to thin disk scale lengths greater than one. Comparing our 
simulation results with observations of the Milky Way and nearby galaxies shows that our results for minor mergers are consistent 
with observations of the ratio of thick to thin disk scale lengths and with the "Toomre diagram" (the sum in quadrature of the vertical 
and radial versus the rotational kinetic energies) of the Milky Way. The simulations of clumpy disks do not show such agreement. 
We conclude that minor mergers are a viable mechanism for the creation of galactic thick disks and investigating stars at several kpc 
above the mid-plane of the Milky Way and other galaxies may provide a quantitative method for studying the (minor) merger history 
of galaxies. 

Key words, galaxies: interaction - galaxies: formation - galaxies: evolution - galaxies: structure and kinematics 
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1. Introduction 

Ever since the existence of a t hick disk component w as shown by 
\ star counts in the MiUcy Way dGilmore & Reidll983h and by stel- 
. lar lumin osity distribution studies of external early- and late-type 

■ galaxies ('Burstein"l979'; 'Tsikoudi''l979'; van der Kaiit & Searld 
il981a; Shaw & G ilmore 1989; Dalcari ton & Bernstein 2002 

■ iPohlen et al .1 120041: IXikhonov et al.li2005l) . its origin and proper- 
' ties have been the subject of a continuing vigorous debate. The 

distinctive properties of their stars, e.g. their metallicity distri- 
bution and kinematics, suggest that the thick disk is a distinct 
and perh aps inte rmediate c omponent between the thin disk and 
the halo (lNorrisl[r9 86: Wvs e & Gilmorelll986llSandage & FoutsI 
ll987HCarnev et al.iil989i:,Ratnatunga & Freemanlll9'89lK 



or from directly accretin g dynamically hot star s throu 



mergers 

'I989;'0uinn et al 



and interactions ( Statler 19881 



Herngui st & OuinnI 
AgueiTi et al. 200^^ 
2005i : ,Brook et all 



A large variety of models has been proposed to ex- 
plain the formation of thick disks - involving either secu- 
lar heating of the disk through scattering by spiral waves 
jBarbanis & Woltierl967:Carlberg & Sell wood,1985j, in olecu- 
lar clouds (ISpitzer & S chwarzschild 195 11: lLacevlfl984 . bars, 
self-gravitating clumps (Boumaud et al., ,2009L or heating disk 
stars through external processes such as galaxy mergers, in- 
falling clumpy gas coalescing into a thin disk, and proto- 
galactic fragments coalescing before thin disk formation. On 
the other hand, a thick disk component can also be formed 
by violent heating and relaxation due to satellite accretion(s) 



1993: Walker etal. 1996: 

l Abadi et al.. 2003: ,Yoachim & Dalcanton 

I2OO7I; IVillalobos & Helm il l2008l l2009irA s recentlv discussed 
bv ISales et al.l (l2009h and iDi Matteo et alJ (1201 Ol) . these differ- 
ent formation mechanisms should produce different signatures in 
the eccentricities of the stellar orbits, thus providing a potential 
diagnostic to disentangle the dominant formation scenarios (see 
iDi Matteo et"ani2010l: iPierickx et al.ll20Tol: IWilson etaDl2010l 
and Casetti-Dinescu et al. for a comparison of model pre- 

dictions and observations). Nevertheless, the main astrophysical 
processes which drive thick disk formation remain unclear. 

The structure, kinematics and enrichment of the thick disk 
provide important clues for solving the mystery of its origin. 
Many photometric and spectroscopic observations show that 
thick disk stars are generally old, and have lower metallic- 
ity and larger v e locity dispersion than thin disk stars issenl 



1995 



2003 



IChiba & Beers"2 ()0d iGilmore et aDl2002l: ISoubiran et al 
Parker et al.. ,20041) . The vertical surface density pro 



file of thick disk stars follows a sech" profile with a scale 
height several times larger than the thin disk, whereas their 
radial profile has a scale length that is generally 1 to 2 
times larger than that of the thin disk ( Morrison et aTl 19971 
lAbe et all [T999t iDalcanton & BemsteinI 120021: iNeeser et al.l 
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l2002HPohlen et alJl2004t lYoachim & Dalcantonll2005h . In mod- 
els, the thick disk scale height is usually assumed to be in- 
depen dent of radius (Pohlen et al. 2004; Yoachim & D alcantonI 
I2005h . in accordance with earlier observational studies on 
thin and thick disk scale heights ( van der Kruit & Searli 



T98Tallbt IShaw & Gilmord \WM iNaslund & Jorsated 119971 



Rauscher et al. However, for a samp le of 48 early- and 



late-type galaxies de Grijs & Peletierj (1 19971) found an increase 



in disk scale height with radius, where early-typ e galaxies show 
the largest increase. Also lNaravan & Jog|(l2002h pointed out that 
an increasing s cale height was consistent with tw o of the galax- 
ies studied by Ivan der Kruit & Searig (Il981albl) . and that their 
vertical velocity dispersion distribution does not scale as twice 
the disk scale length with radius, thus supportin g the possibility 
of a m oderate flaring in these disks. Recently, iBournaud et al.l 
(I2009h modeled the formation of a thick disk through gravita- 
tional instabilities and pointed out that thick disks produced in 
this way have constant scale heights, while this is not likely to 
be the case if they were formed through minor mergers. 

Therefore, the properties of vertical disk profiles can 
help in understanding their formation processes. It ap- 
pears that t hick disk scale heig ht is related to both 

alaxy mass (lYoachim & Dalcanto nr i2006l) and Hubble type 



Moreover, some disk-dominated 
alaxies do not show an obvious thick disk, e.g. NGC 4244 



Ide Griis & Peletiei 



lFrvetal.1 [19991 but seelc^omeron et"dll2011l for evidence of 
a subtle thick disk), and some even have thick disks that 
counter-rotate relative to the ir thin disks ([Morrison et al 
lYoachim & Dalcantonll2005l) . These observations, taken all to- 
gether, fav or a merger- induced thick disk formation mechanism 
(see, e.g., lYoachim & Dalcanton 2005, 200(i and references 
therein). 

Int erest ingly, it has r ecently been shown by iPurcell et al.l 
(1201 Ol) and IZolotov et al.l ([2010) that minor and major mergers 
cannot only efficiently heat a pre-existing stellar disk, but also 
contribute to the building of galaxy halos. A trace of this heated 
population may be found in the inner Milky Way halo, as re- 
cently discussed bv lNissen & Schusteii (120101) . 

Despite considerable effbrt, we are still far from having a 
detailed understanding of the origin and properties of galactic 
thick disks. In this paper, we investigate if N-body/SPH numeri- 
cal simulations of minor mergers, with a mass ratio of 10:1 and 
20: 1 , can produce realistic thick disks. We concentrate on certain 
aspects of post-merger thick disks - i.e., vertical surface density 
profiles and disk scale heights - that are well suited for compari- 
son with the observed properties of thick disks in the Milky Way 
and nearby galaxies. In particular, we investigate the "excess" of 
stars at great heights (z>2 kpc) from the galaxy mid-plane that 
is naturally formed thro ugh minor mergers. This stellar excess, 
as alre ady discussed by iPurcell et al.l (1201 Ol) and IZolotov et al.l 
([20Toh . may have contributed to the stellar population currently 
found in inner galactic halos. For the first time, we describe its 
main properties and the impact on them of initial orbital con- 
figurations, primary galaxy gas fraction and consecutive minor 
mergers. In addition, we compare the scale length ratios of the 
thick and the thin disk to observations, and the Toomre diagrams 
of minor mergers with that of the Milky Way. Furthermore, we 
will compare the vertical structure of thick disks formed through 
minor mergers and through scattering by clumps in the distri- 
bution of mass within the thin disk, and discuss the expected 
differences in their observable properties. 

The paper is organized as follows. In §[2] we describe the 
numerical code, the initial galaxy models and orbital conditions 
adopted for simulation runs. § [3] presents the vertical structure 



and kinematics of merger-induced thick disks and their depen- 
dence on initial configuration parameters in our simulations. 
Finally, we discuss the results in § |4| and draw our conclusions 
in §[3 



2. Models and initial conditions 

We study the interaction and coalescence of a satellite galaxy 
with a much more massive disk galaxy. These minor merger 
simulations are part of an on-going program to model and un- 
derstand the role minor mergers play in the evolution of an- 
gular momentum an d morphologica l and kinematic properties 
of galaxies (see also lOu et al ] l2010l) . The massive disk galaxy 
(primary galaxy) consists of a dark matter halo which is ini- 
tially not rotating, a central bulge, and a stellar disk - hereafter 
we refer to this model with the nomenclature gSO (=giant SO 
galaxy). We also consider initial models containing a gaseous 
disk, whose mass, Mg„s, is 10% or 20% of that of the stellar disk, 
and to which we refer as "gSa" and "gSb", respectively. The 
spherical dark matter halo and the stellar bulge are represented 
by Plummer spheres (iBinnev & Tremainell 19871) and neither has 
any rotation at the beginning of the simulation. The stellar and 
gase ous disks are represented b y Miyamoto-Nagai density pro- 
files (IBinnev & Tremainell987l) . Table [T]lists the total masses of 
the halo, bulge and disk components (Mh, Mb, M, and Mg„s) 
and the core radii of the halo and the bulge (rn and r^), as well 
as the vertical and radial exponential scale lengths of the stellar 
and gaseous disks (h,, a» and hgas, Ogas)- The elliptical satellite 
galaxy has a total mass equal to 10% or 5% of that of the primary 
galaxy - we refer to these models as dEO and sEO, respectively. It 
consists of a spherical stellar and a dark matter component, both 
modeled with initially non-rotating Plummer profiles, whose pa- 
ramete rs are also given in Table [1] As we have done in lOu et al.l 
(1201 ol) . all galaxy models are evolved in isolation for 1 Gyr be- 
fore the interaction starts. The orbital parameters for these mi- 
nor interacti ons have been de scribed in lChilingarian et all (1201 Ol 
Table 9) and lOu et al.1 (1201 Oi Table 2) and we refer the reader to 
these papers for a detailed description. 

In th i s pap er, we follow the nomenclature adopted in 
lOu et al.l (I2OIOI) . i.e., using a six-character string to indicate 
the morphology of the interacting galaxies: the first three de- 
scribe the type of the primary galaxy, "gSO", "gSa" or "gSb", 
depending on its gas fraction, and the following three describe 
the satellite galaxy, dEO or sEO. This is followed by the suf- 
fix "dir" or "ret", forprograde or retrograde orbits, respectively. 
The subsequent two numbers, either "33" or "60", indicate the 
initial inclination of the satellite orbit with respect to the disk 
plane of the primary galaxy. Thus, for example, the nomencla- 
ture "gS0dE001dir3 3 "refers to a prograde encounter between a 
gSO galaxy and a dwarf elliptical satellite galaxy, with an ini- 
tial inclination ; - 33°, whose initial orbital parameters are 
those co rresponding to th e id="01dir" case listed in Table 9 of 
Chilingarian et al.l (I2010h . 

To study the efi'ect of multiple mergers on the properties of 
post-merger stellar disk, we ran some simulations in which the 
primary galaxy accretes consecutively two satellites over a pe- 
riod of 3-5 Gyr We have considered both the case in which both 
satellites coming from a fixed direction inclined 33° with respect 
to the primary disk, as well as the case in which the second satel- 
lite comes from a different direction - hereafter we will refer to 
these models as "multi-A" and "multi-B" respectively. For exam- 
ple, the nomenclature "multi-A, 2x10: 1" refers to the successive 
mergers of two satellite galaxies whose mass is one-tenth of that 
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Table 1. Parameters of the initial models of halo, bulge, stellar 
and gaseous disks. 





gSO 


gSa 


gSb 


dEO 


sEO 


Mb [2.3 X 10' Mo] 


10. 


10. 


5. 


7. 


3.5 


M« [2.3 X 10' Mq] 


50. 


50. 


75. 


3. 


1.5 


M. [2.3 X 10' Mo] 


40. 


40. 


20. 










0.1 


0.2 






tb [kpc] 


2. 


2. 


1. 


1.3 


0.9 


th [kpc] 


10. 


10. 


12. 


2.2 


1.55 


a, [kpc] 


4. 


4. 


5. 






h, [kpc] 


0.5 


0.5 


0.5 






flj„j [kpc] 




5. 


6. 






hga, i:kpc] 




0.2 


0.2 







Table 2. Particle numbers for primary galaxies and satellites 





gso 


gSa 


gSb 


dEO 


sEO 


N 

gas 
^ star 

Ndm 


320,000 
160,000 


80,000 
240,000 
160,000 


160,000 
160,000 
160,000 


32,000 
16,000 


16,000 
8,000 



of the primary galaxy and whose orbits are inclined 33° with 
respect to the primary disk. 

All simulatio ns (50 in total) were run u sing the Tree-SPH 
code described in lSemelin & Combesl(l2002l) . A total of Ntot - 
528, 000 particles have been distributed between the primary and 
satellite galaxies (see Table |2]i for all simulations. We also test 
the dependence of results on the particle number used in the sim- 
ulation by running some additional simulations with a total of 
Ntot = 1,056,000 particles. A Plummer potential is used to 
soften gravity on small scales, with constant softening lengths 
of e = 200 pc (e = 170 pc for the high resolution simulations), 
for all particles. The equations of motion are integrated using a 
leapfrog algorithm with a fixed time step of 0.5 Myr. With these 
choices, the relative error in the conservation of the total energy 
is about 10"^ per time step. 

In this study we also compare the vertical structure of mi- 
nor merger-induced thick disks with those generated by an inter- 
nal process, through scattering by massive clum ps formed in an 
initiall y unstable disk, as recently proposed bv iBoumaud et al] 
(I2009h . For this, we anal yzed two simulat i ons of gas-rich, un- 
stable disk galaxies from iDi IVIatteo et al.l (l2008h . These use a 
total number of 120,000 particles, equally distributed among the 
gas, stars, and dark matter components, whose parameters are 
the same as those of the gSb model (see Table [1]), except for the 
gas mass, which is initially 50% that of the stellar disk. These 
gas-rich galaxies are characterized by a very small initial Toomre 
disk stability parameter for the gaseous component [Qgas = 0.3 
and Qgns - 0.1, hereafter called "gSb-Hul" and "gSbH-u2", re- 
spectively). Reducing the local disk stability causes the gaseous 
disk to fragment into many clumps, particularly in the "gSb+u2" 
case. Once formed, these clumps migrate towards the galaxy 
center; some dissolve under the influence of the tidal field of the 
disk, while others survive and merge in the galaxy center, con- 
tributing to the formation of the central bulge. As recently pro- 
posed by Bournaud et al. (2009), these clumps can scatter stars 
in the pre-existing thin disk, contributing to the formation of a 
thick disk component. 



3. Results 

During the merging process, tidal interactions between merging 
galaxies deposit part of their energy into stellar motions in the 
disk of the primary galaxy, thus causing the disk to become kine- 
matically hotter and to spread out both radially and vertically, 
changing its morphology as well as its kinematics. In this sec- 
tion we will analyze both dissipationless and dissipative minor 
mergers and their possible role in the formation of the thick disk. 

3.1. Characteristics of tlie post-merger stellar disk 

Here we investigate the vertical stellar distribution of remnant 
disks in minor mergers, for which we assume a generalized sech 
profile, X(r, z) = I.rSech-^"{--^), where is the radial expo- 
nential surface density profile and zo is the disk scale height. 
This profile can fit, for example, an isothermal sheet if n = 1, 
and an exponential if n - oo - the two limiting cases of this 
generalized profile. Base d on the study of the three-dimensional 
distribution of disk stars, Ivan derKruiddlgsl found that an in- 
ter mediate disk model with n - 2 is the most appropriate (see 
also [Schwarzkopf & De ttmar 2000). whi ch is a lso supported by 
the theoretical study of iBaneriee & Joel (l2007b . We analyze, as 
far as possible, the properties of remnant disks as they typically 
would be observed. First, we project the post-merger stellar disk 
to edge-on, then average the vertical surface density profiles on 
both sides of the stellar disk mid-plane and fit them using a sech 
profile with « = 2, i.e., sech{-^), allowing the scale height pa- 
rameter Zo to be unconstrained. We also fit the disk scale height 
for the two limiting cases of the generalized sech profile, n - oo 
and « = 1: the exponential fit gives an on average 1.4 times 
larger disk scale height than the sech profile, whereas the sech^ 
fit gives a value which is 0.8-0.9 times smaller. 

The vertical stellar surface density profiles are analyzed be- 
tween 0.5 and 1 Gyr after the merger is complete, at a time when 
the post-merger stellar disk has settled roughly into a new equi- 
librium configuration. At later post-merger times, i.e., f > 1 Gyr, 
the difference in disk scale height with respect to the value at 
t - 0.5 - 1 Gyr is no more than 15% inside r < 6r,;. The stellar 
disk scale length, r^, has been evaluated by projecting the galaxy 
to face-on. The diff'erence between the final and initial disk scale 
length is always small, Arj < 10%, except for dissipative mi- 
nor mergers with a gas-to-stellar mass fraction /go, = 0.2, where 
the disk scale length shrinks by ~ 15 - 20%. Thus the exact 
time at which the merger simulations are analyzed beyond 0.5- 
1 Gyr after the completion of the merger makes little difference 
in the final results. In this study we also mask the contribution of 
the bulge component by excluding stars that were initially in the 
bulge of the primary galaxy and we consider only stars which 
were initially in the disk of the primary galaxy or in the satellite. 

We note that in the case of dissipative minor merger simula- 
tions, a new stellar disk forms during the merging process due to 
star formation in the gaseous disk. This "new" stellar disk com- 
ponent is thin, with scale heights of 300 - 400 pc only (see Fig. [T] 
as an example). We will not discuss in detail this newly formed 
thin disk component, but concentrate on the old stellar disk, i.e., 
the one that was akeady in place before the interaction, in order 
to make a direct comparison with that found in dissipationless 
mergers. We do not find any young stars at scale heights larger 
than that of the thin disk formed during the merger. 

In the rest of this section, we will examine the vertical struc- 
ture of the remnant disk after a single dissipationless minor 
merger, to illustrate the process of thick disk formation during a 
gas-free minor merger. We will present the vertical stellar distri- 



initial primary stars 
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Fig. 1. Stellar distribution in a gSb galaxy model in its initial configuration (left panel) and at 3 Gyr after the start of the simulation 
(middle and right panels). All panels show the galaxy from an edge-on perspective. Of the post-merger galaxy, both the original 
stellar population and the new thin disk formed during the merger are shown in the middle panel, while in the right panel we 
show only stars from the original stellar population of both galaxies. The thickening of the disk is entirely due to stars from the 
original thin disk and the satellite. The same color and contour scales are used in all panels. Surface density contours are spaced 
logarithmically and each panel is 40 kpcx40 kpc in size. 



bution and the fitting results, and will later compare these prop- 
erties with those of remnant thick disks found in dissipative or 
two consecutive minor merger simulations, as well as with thick 
disks formed through massive clumps in a gravitationally unsta- 
ble disk. 

3.1 .1 . Vertical structure of the post-merger thick stellar disk: 
evidence for a stellar excess 

The vertical surface density profile of a thick disk formed 
through a dissipationless minor merger is shown in Fig. |2] This 
profile has been evaluated 0.5 Gyr after the completion of the 
merger The main features of the profile fits are: 

- A single sech function cannot fit the entire vertical profile. 
While the lower regions, at z<2 kpc, can be described by 
a sech function, after removing this sech function there re- 
mains a large excess of stars at greater heights above the 
mid-plane (z>2 kpc); 

- This excess is especially prominent in the inner regions of 
the disk, at r < 3r^, and it virtually disappears at larger radii; 

- The excess has a scale height several times larger than that 
of the main thick disk component (i.e., the one at z<2 kpc): 
at r = 0, for example, the scale height of the stellar excess 
is zo ~ 2.3 kpc, about 5 times that of the main thick disk 
component, which has zo ~ 0.42 kpc; 

- The excess has a scale height that is approximately con- 
stant with distance from the galaxy center, while the thick 
disk scale heigh t incre ases with radius, as already sho wn by 
iBoumaud etall (l2009l) and iKazantzidis et alJ (l2008h - this 
will be discussed further in § 13.21 

- Due to this excess, two sech functions of different scale 
heights are needed to fit the entire vertical stellar profile. 



3.1 .2. Thick disk and stellar excess: formation, growth and 
composition 

How do the thick disk and the stellar excess develop in the sim- 
ulations? We show an example of the formation and growth of 
a thick disk in Fig. [3] After the first pericenter passage of the 
satellite galaxy ( at f = 0.5 Gyr) the thin stellar disk compo- 
nent has already thickened and grown vertically, and after the 
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Fig. 2. Vertical stellar surface density profiles of a post-merger 
thick stellar disk at 0.5 Gyr after the completion of a dissipa- 
tionless minor merger (black). We show profiles at a variety of 
distances from the galaxy center along the galactic mid-plane, 
in units of the disk scale length rj. Also shown are the double 
sech functions fitted to the profiles and their corresponding scale 
heights Zo in kpc (labeled in red). 



second pericenter passage, at f = 1 .5 Gyr, the stellar disk has 
acquired the double sech profile. Since the stellar excess is al- 
ready in place well before the satellite merges with the primary 
it is not associated generally with stars from the satellite galaxy, 
but rather composed of stars which were initially in the primary 
disk and subsequently heated, reaching greater heights during 
the merging process. In all dissipative mergers and in most of 
dissipationless ones, stars initially in the di sk dominate at all disk 
heights, in agreement with the results of iPurcell et al.l (l2010l) . 
Only in few dissipationless mergers do we find that the outer 
vertical profile is dominated by stars from the satellite galaxy. 

The formation of a double sech profile in the vertical density 
distribution is not due to the presence of the bulge. Although we 
excluded bulge stars from the analysis of the vertical profile of 
the thick disk, we did run simulations of dissipationless minor 
mergers with bulgeless galaxies. In these simulations we also 
found a stellar excess with characteristics very similar to simu- 
lations including a bulge component. 
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Fig. 3. Upper panel: Vertical stellar surface density profiles as 
function of radius and time (indicated by color; see lower panel) 
for a dissipationless minor merger The radius is indicated by the 
4 black labels. Lower panel: The distance between the centers of 
the satellite and the primary galaxy as function of time. Arrows 
mark the merger phases which correspond to profiles with the 
same color and line styles shown in the upper panel. 



In summary, stars in the stellar excess have their origin in dif- 
ferent regions of the thin disk, which means that they are not only 
heated vertically, but also redistributed radially. As discussed in 
§ m this has consequences for the expected metallicity gradient 
in the merger remnant. 



scale heights increased with disk radius within 4 or 5 disk scale 
lengths, and that the hvei/fd ratio ranged from about 2 to 4. 

As pointed out by Ouinn et al. (1993), in mergers of pri- 
mary galaxies with their satellites, the stellar disk is heated in 
three dimensions and although the velocity dispersion profiles 
of the remnant disk are similar to those of the initial disk, they 
show more extended exponential profiles. Fitting both the initial 
and post-merger vertical velocity dispersion profiles with expo- 
nential functions, we found that while the initial profiles usu- 
ally have hy^i/rj ~ 2, in agreement with an initial stellar disk 
of constant scale height, as assumed in our model, the hyei/fd 
ratio of the post-merger vertical velocity profile was increased 
to the range of ~ 2. 5 - 4, which is in agreement with what 
iNaravan & JoS ( |2002|) found. 

From the definition of zo deduced from Poisson's and Jeans' 
equations, this link between the hy^i/rd ratio and the thick disk 
scale height is not difficult to understand. Assuming the velocity 
dispersio n follows an exponential profile , the definition of zo (see 
Eq.(3) in Ivan der Kruit & Sear i3ll981bl) can be rewritten as 



Co exp 



hyd / 2rd 



(1) 



where Co = 



o-;('-=0) 



„,, 7T and G is the gravitation constant. 
cr,{r - 0) and po are the vertical velocity dispersion and the 
volume density at the galactic center, respectively. It is quite 
clear from Eq.(l) that zo is independent of radius if hyei/rd = 2, 
whereas zo increases with radius when hyei/rd > 2. 

Fig.|4]also shows the scale height of the stellar excess present 
in the vertical surface density profiles seen at z > 4 - 5 kpc and 
^ 3rd (see also Fig. |2]i. The final scale height of the stellar 
excesses in the simulations are independent of radius and are 
sensitive to the characteristics of the satellite orbit. Variations in 
the orbital parameters lead to a large scatter in the vertical scale 
height of the stellar excesses, of up to 30% i\Azo\/ zo,mean ^ 30%). 



3.2. Varying orbital and morpliological parameters 

3.2.1. Dependence on initial orbital configurations 

Regardless of the initial orbital parameters (orbital energy and 
angular momentum, direct and retrograde orbits, initial orbital 
inclination of the satellite orbit) we always find a thickened stel- 
lar disk in the merger remnant, whose scale height increases with 
the radial distance along the major axis of the disk-plane (Fig.|4l). 
For the dissipationless mergers over the range r^ < r < 3r^/ the 
scale height is zo = 1-4 - 2.5 kpc, or about 3-6 times larger than 
the Zo value near the disk center, where zo = 0.4 kpc. The six 
sets of diff'erent orbital configurations produce a scatter around 
the mean value of |Azol/zo,Hiean ^ 15%. The increase in the disk 
scale height with radius implies that the disk heating is the most 
significant in the outer regions of the disk. We do not find such 
a strong thickening of the stellar disk if the galaxy is evolved in 
isolation for 3 Gyr (see Fig. |4]i: the disk scale height increases 
only marginally, less than 25% at r ~ rj, and the thickening is 
independent of radius. 

Early research in self-gravitating stellar disks suggested that 
the thick disk scale height is independent of radius, assum- 
ing that the vertical velocity dispersion decreases exponentially 
with radius and relates to the disk scale length as hy ei - 2rd 
(Ivan der Kruit & Searlelll981ah . lNaravan & Jogl(l2002l) general- 
ized this picture by taking into account the gravitational contri- 
butions from the gas and the dark matter halo and allowing the 
ratio hyeilrj to vary. Contrary to previous results, they found that 



3.2.2. Dependence on satellite central density 

Do the flaring of the thick disk and the scale height of the stel- 
lar excess found in Fig. |4] depend on the central density of the 
satellite galaxy? To investigate this point, we ran some simula- 
tions changing the central density of the baryonic component of 
the dEO satellite, using satellites having a central volume den- 
sity 50% higher (dEOh) or lower (dEOl) than the reference dEO 
galaxy. A denser satellite heats the thick disk more effectively, 
especially in the inner regions of the primary galaxy (inside 
~ 2r,;), thus resulting in a thick disk scale height about 20% 
higher than that found after the merger of less dense satellites. 
However, no clear dependence on the satellite central density is 
found for the characteristics of the stellar excess (see Fig.|5]l. 

3.2.3. Dependence on gas fraction 

How does including gas in the primary galaxy's disk affect the 
outcome of these simulations? In Fig. |6] we compare the disk 
scale heights of dissipative minor mergers (with gas-to-stellar 
mass fractions fgas =0.1 and 0.2) with those of dissipation- 
less mergers. A thick stellar disk is present in all merger rem- 
nants, independent of the gas fraction in the primary, and its 
scale height always increases with radius. Dissipative mergers 
with fgas =0.1 have similar scale heights as the dissipationless 
(gas-free) ones. However, in mergers with fgas - 0.2 a signifi- 
cantly thinner remnant stellar thick disk is formed, of which the 
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Fig. 4. The scale heights zq of the merger-induced thick disk and 
stellar excess as function of radius in units of disk scale length, 
r/rd, for dissipationless minor mergers with six sets of different 
initial orbital parameters (see also FigjS). The shaded regions in- 
dicate the scatter in the scale heights. Most of this scatter results 
from the differences in initial orbital configurations. Also shown 
are the scale heights of the original stellar disk (dotted line) and 
after its evolution in isolation for 3 Gyr (dashed line). 
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Fig. 5. The scale heights of the merger-induced thick disk and 
stellar excess as function of radial distance in units of disk scale 
length, r/r^, for dissipationless minor mergers with three differ- 
ent baryonic central mass concentrations of the merging satellite. 
The mass of the satellite is the same, but the central volume den- 
sity is 50% higher (dEOm, dashed lines) or 50% lower (dEOI, 
dotted lines) than in the reference satellite dEO (solid lines). 



scale height at r > 2rj is on average ~20%-30% smaller than in 
dissipationless mergers. Interestingly, this reduction is compara- 
ble to the 2 5% found for 10:1 d issipative mergers with 20% gas 
fraction by iMoster et al.l (1201 Ol) . which indicates that their con- 
clusions are robust to the use of different initial galaxy models, 
orbital conditions and numerical codes. 

On the other hand, the stellar excess in dissipative minor 
mergers still has scale heights that are independent of radius, and 
very similar to those found in dissipationless mergers (see left 
panel in Fig.|6]l. This implies that the dissipative gas component 
and further induced star formation mainly affects the higher den- 
sity regions close to the remnant disk mid-plane and that it has 
almost no effect in regions at greater heights above the disk mid- 
plane. This is consistent with the excess being predominately 
due to tidal effects early in the merger event. 




stellar excess 



2 
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Fig. 6. Scale heights of the merger-induced thick disk and stellar 
excess as function of radial distance in units of disk scale length, 
r/r,;, for minor mergers with three different gas-to-stellar mass 
fractions, fgas = 0, 0.1 and 0.2. The shaded regions indicate the 
scatter in the scale heights due to the variety of initial orbital 
configurations. 



3.3. Dissipative minor mergers versus clump instabilities 

It is of great interest to investigate the differences in the prop- 
erties of merger-induced thick disks to those formed by secular 
processes such as scattering due to clumps for med in a gravita- 
tionall y unstable disk, as recently suggested by iBoumaud et al.l 
(l2009l) . To this end, we compare the vertical surface density pro- 
files of thick disks induced by minor mergers with those formed 
in gas-rich galaxies which are initially unstable to clump for- 
mation and were evolved in isolation for 3 Gyr (see §|2]for de- 
tails of these galaxy models). We found that the resulting scale 
height in the clumpy isolated galaxies is about 2 times larger, 
Zo ~ 0.85, than the initial val ue and nearly constant with radius, 
thus confirming the results of iBournaud et al.l (l2009h . Scattering 
by gas clumps do not only produce thick disks with constant 
scale heights, but also vertical surface density profiles which do 
not show any stellar excess (Fig. |7]i. In principle, this may indi- 
cate an observational way to distinguish a secular from a merger 
origin for thick disks. We will elaborate on this further in the §|4] 



3.4. Consecutive minor mergers and their impact on the thick 
disk and stellar excess scale heights 

In hierarchical structure formation scenarios, galaxies grow 
through both minor and major mergers and periods of gas ac- 
cretion from either cooling gas in their halos or from cosmic fil- 
aments. Satellite accretion is expected to be quite common in a 
ACD M universe and important for galaxy growth (iGuo & White! 
l2008l) . The dynamical timescale for a satellite in the dark mat- 
ter halo of a significantly more massive primary galaxy can be, 
and often is, long compared to the merging time scale. In such 
a situation we may expect multiple minor mergers to be com- 
mon since this allows the primary galaxy to have several satellite 
galaxies that have yet to merge. 

What is the effect of consecutive mergers on the vertical stel- 
lar distribution of disks and their impact on the thick disk and 
stellar excess? Do the effects "saturate" - in other words, does 
the additional heating caused by each subsequent merger have 
an increasingly smaller effect on heating of the disk? In order to 
gain insight in this, we compared the impact on the heating of 
the disk of a single minor merger to that of two sequential minor 
mergers. We adopted four merger histories and two mass-ratios 
- a single 20: 1 merger and two sequential 20: 1 mergers, a single 
10:1 and two sequential 10:1 mergers, which we will label as 
20:1, 2x20:1, 10:1, and 2x10:1, respectively. All of these merg- 
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Fig. 7. Vertical density profiles of the post-merger stellar disk 
in a dissipative minor merger with fgas - 0.2 (upper panel) 
and of an isolated "clumpy disk" galaxy with fgas - 0.5 (lower 
panel). Profiles are shown at four different exponential disk scale 
lengths, r^i, from the galaxy center along the major axis of the 
galactic mid-plane. The best fitting sech(z) functions to these 
profiles and corresponding scale heights are also shown (in red). 
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Fig. 8. Vertical stellar surface density profiles of a post-merger 
thick stellar disk 0.5 Gyr after the completion of a single dissi- 
pationless minor merger with a mass ratio of 20: 1 (upper panel), 
and two consecutive dissipationless minor mergers with a mass 
ratio of 10:1 (lower panel). We show profiles at four different 
exponential disk scale lengths, rj, from the galaxy center along 
the major axis of the galactic mid-plane. Superposed on these 
distributions are the fitted double sech functions and their corre- 
sponding scale heights (in red). 



ers are dissipationless and intended only as an initial exploration 
of the effects of multiple minor mergers on galaxy disks. 

The vertical surface density profiles of the remnant of a 20: 1 
and a 2x10:1 merger are shown in Fig. |8] Comparing these pro- 
files with those of a single 10:1 merger (Fig.|2}, we find, as ex- 
pected, that the thickness of the final disk increases with the total 
mass acquired by a galaxy through mergers. In other words, the 
more mass accreted, the more internal energy accrued, and the 
more efficient the heating of the stellar disk. But changing the 
mass of the satellite has an impact also on the formation and re- 
sultant scale height of the stellar excess. A 20:1 minor merger 
bring less energy to the system, which means not only a smaller 
thick disk scale height, but also a obscure stellar excess: stars 
originally in the thin disk are not heated enough to reach high 
vertical distances from the galaxy mid-plane. On the other hand, 
a 2x 10: 1 heats the thick disk considerably, especially in the outer 
parts. This depends on our original definition of thick disk and 
stellar excess. The behavior of thick disk and stellar excess scale 
heights as a function of radius, for six representative single and 
two consecutive mergers, is shown in Fig.|9] 

In the outer regions, beyond Irj, the scale height of the rem- 
nant disk in the single 10:1 minor merger case is ~1.5 times 
larger than in the 20:1 merger, whereas the 2x20:1 merger 
has a similar disk scale height as a single 10:1 minor merger 
Varying the orbital inclination of the second satellite ("multi- 
A" and "multi-B" cases) only causes a rather trivial change in 
the disk scale height. The effect on disk scale height of two 



20:1 mergers appears to be cumulative, and almost the same 
as for a single 10:1 merger (Fig.|9). Thus comparing a 2x10:1 
merger is like comparing with a 4x20: 1 merger, provided the ef- 
fect of multiple mergers does not saturate. From the very simple 
arguments of Liouville's theorem and the conservation of en- 
ergy, this is expected to be the case for dissipationless systems 
(Binney & Tremaine 1 987) . 

In our simulations of a 2x10:1 minor merger, after the first 
10:1 merger the resultant thick disk still responds dynamically 
to the heating induced by the subsequent 10:1 merger (Fig.|9). 
After the second merger the final disk sc ale height is 1.5-2 times 
larger than in a single 10: 1 merger Similarly, the stellar excesses 
found at greater heights in the vertical stellar distribution have 
also increased in the subsequent merger events (Fig. |9) - their 
disk scale heights after a single 10:1 merger are 1.5 times larger 
than that after a single 20:1 merger In the 2x10:1 mergers, their 
scale heights are not only ~ 1 .4 - 2 times larger than in single 
10:1 mergers but they also increase with radius. Thus, we find 
that the effect of two consecutive minor mergers in dissipation- 
less systems does not saturate but is cumulative. However, we 
note that after the second 1:10 merger the stellar excess begins 
to develop a scale height that depends (weakly) on radius, and is 
more pronounced at larger radii where its scale height becomes 
very similar to that of the thick disk itself (Fig.|9]l. Thus, at least 
for the highest mass ratios and in the outer disk, the stellar ex- 
cess becomes almost indistinguishable from the thick disk. For 
most of our two consecutive interaction models we confirm what 
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Fig. 9. Scale heights of the merger-induced thick disk and stel- 
lar excess as function of radial distance in units of disk scale 
length, r/r^/, after a single or two consecutive minor mergers. 
We show single mergers with a mass ratio of 10:1 and 20:1, and 
2x10:1 and 2x20:1 consecutive merger events. In the consecu- 
tive merger cases the satellites are either accreted from the same 
direction (multi-A) or from different directions (multi-B). 



we found for the single mergers, that is stars from the primary 
disk galaxy dominate the vertical stellar distribution everywhere 
in the thick disk. 



3.5. Kinematics of tliick disk and stellar excess 

Do the thick disk and stellar excess have similar kinematic prop- 
erties? A detailed study of various kinemati c properties of m inor 
merger remnants will be made separately (lOu et al.ll201 lb . and 
here we will discuss only one of the kinematical characteristics 
of the stellar excess: its rotational lag with respect to thick disk 
stars. This rotational lag is clear in Fig. [10] where the tangential 
velocity of stars is shown as a function of the vertical distance 
from the galaxy mid-plane, for a dissipationless and a dissipative 
merger In both cases the tangential velocity decreases with dis- 
tance from the mid-plane, thus indicating that stars far from the 
mid-plane plane, which populate the stellar excess, have lower 
specific angular momentum than stars in the thick disk (which 
are close r to the mid-plane). Th is is consistent with the trend 
found bv iDi Matteo et all (1201 Ol) for the eccentricities of stellar 
orbits after minor mergers: stars tend to move on more eccentric 
orbits as their vertical distance increases and thus tend to rotate 
more slowly than they did originally. 

This increase in the rotational lag is qualitatively in agree- 
ment with observations (Ivezic et al. 2008), and based on our 
models we would expect to find high metallicity, small lag stars 
at low height, and stars with lower metallicity and larger lag at 
greater height - this is due to the way the thick disk is formed in 
our model, through the heating of a thin disk that continuously 
forms stars. 



4. Discussion 

To what extent do the properties of our modeled thick disks 
formed through minor rnergers agree with obs ervations? As 
previo usly discussed by Purcell et al. (20 13) and IZolotov et alj 
(120 lOl) . minor mergers are able to considerably heat a pre- 
existing thin disk. We have shown (§ 13. Il l that the resulting stellar 
vertical profile is characterized by a thick disk component and a 
stellar excess, both of which can be fitted by sech profiles with 
different scale heights. The excess is a natural characteristic of 
minor merger models, which cannot be reproduced by secular 
mechanisms, not even through scattering by clumps in unsta- 
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Fig. 10. Tangential velocity profiles of the post-merger stellar 
disk as function of height from the galaxy mid-plane in dissi- 
pationless (upper panels) and dissipational (lower panels) minor 
mergers with a gas mass fraction of 0.2. Profiles are shown at 
different disk scale lengths, r^, from the galaxy center along the 
major axis of the galactic mid-plane. 



ble gas-rich disks. Identifying these excess stars and determin- 
ing their properties may thus provide a way to constrain thick 
disk formation scenarios, but they may be confused with stars 
in the inner halo and be difficult to find and study in practice 
dStanwav et alJl2008l:ICarollo et al.ll2010l) . 

Interestingl y, some of the excess st ars may already have 
been observed. iNissen & Schuster! (1201 Ol) reported the presence 
in the solar neighborhood of stars with halo kinematics, but 
a/Fe abundances similar to those of thick disk stars. These stars 
may well be part of the stellar excess component, as we have 
shown (§ 13.51 ) that excess stars lag rotationally with respect to 
thick disk stars. If the Milky Way thick disk formed through 
the heating of a pre-existing thin component by minor mergers 
at early epochs, we should expect to find an inner halo popula- 
tion (the excess), with abundances similar to those of stars in the 
thick disk, but lagging with respect to them. Compared to stars 
originally in the thin disk, stars accreted from satellite galaxies 
may have abundances which are lower (though probably higher 
than those of dwarf galaxies currently orbiting in the Milky Way 
halo) and probably als o a lower angular momentum content (see 
IVillaloboset"ani201flh . In summary, in the merger scenario the 
inner halo stars would probably have a population with abun- 
dances similar to those of the thick disk, lagging with respect to 
thick disk stars, but rotating more rapidly than accreted satellites. 

The Toomre diagram shows the relationship between the to- 
tal radial and vertical kinetic energy of stars and the rotational 
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Fig. 11. Toomre diagrams of stars at 2-3rd after a prograde minor merger (left), a retrograde minor merger (middle) and in a clumpy 
disk evolved in isolation for 3 Gyr Plotted are the quadrature sum of vertical and radial velocities (u and w) as function of the 
rotational velocity, v, relative to the Local Standard of Rest. The primary galaxy in the merger models has a gas-to-stellar mass ratio 
of fgas = 0.2. The dashed line m arks a total velocity of 180 km s which serves as a potential criteria to separate thick disk stars 
and halo stars (IVenn et al.ll2004h . and the dotted line is the Local Standard of Rest in these models. Also shown are the observed 
thi ck disk stars (triangles), ha lo stars with high a/Fe abundance (squares) and low a/Fe value (crosses) in the solar neighborhood, 
see lNissen & Schusteil(l2010h for further details. The halo stars in retrograde motion are shown in green whereas those in prograde 
motion are in blue. 



energ y relative to the Local Standard of Rest (ISandage & FoutsI 
Il987h . As it is based on kinematic energy, it provides us with an 
appropriate tool to evaluate the results of our models, which are 
kinematic in nature, and a way to compare with the kinematic 
properties of various components in the Milky Way. 

In minor merger models the stellar excess stars lie exactly 
in the area of the Toomre diagram where a enhanced stars are 
found (|Nissen & Schuster 2010), thus strengthening the hypoth- 
esis that they may originally have been in the thin disk. Of 
course, the dynamical and chemical properties of thick disks and 
halos will depend on the properties of the satellite galaxy. We ex- 
pect the scatter in the dynamics and metallicity to increase with 
vertical distance, as the fraction of satellite stars increases with 
height. On the other hand, the clumpy disk model does not re- 
produce the Toomre diagram of the Milky Way (see the right 
panel of Fig.fTTTi. since it cannot heat disk stars as efficiently in 
the vertical direction as minor mergers do. It should be pointed 
out that our models agree with the Toomre diagram of the Milky 
Way, although they are neither tuned to reproduce the Milky Way 
galaxy, nor meant to capture the complex accretion history that 
it may have had. 

Adding a dissipative component (gas) decreases both the 
vertical scale height and the radial scale length of the thick 
disk created during a minor merger (see § 13.3b . For a gas-to- 
stellar mass fraction fgas - 0.2, the decrease in thick disk scale 
height is about ~20% with respect to the dissipationless case, 
and the decrease in scale length is up to 15%-20%. The ratio 
of thick disk scale length and particle-number-averaged scale 
height {rjlzQ,aver) is shown in FiglT2]as function of the circular 
velocity of the disk, Vf,>, for simulations with a variety of mass 
ratio, gas fraction and merger histories. Here iv /r is estimated at 4 
times the disk scale length, r,/, where the disk rotation curve be- 
comes flat, and zo.mcr is the number-of -particles-weighted aver- 
age thick disk scale height within r < A-rj. Since all our primary 
galaxies have the same mass, as we are basically concerned with 
primary:satellite mass ratios in our merger simulations, we only 
sample a single circular velocity. However, we note that more or 
less massive primaries would produce similar results. We com- 
pare our results with observations of late-type galaxies and th e 
MW dLarsen & Humphrevsl2003l: lYoachim & Dalcantonll2005l) . 
For all our minor merger models, whether gas-free or with gas , 



or single or double merger events, the ratios are consistent with 
those of extragalactic thick disks and the Milky Way old thick 
disk. Without a gas component, two consecutive mergers ef- 
fectively halve the ratio by doubling the disk scale height and 
keeping the radial scale length almost unchanged. The impact 
of accreting additional gas on the scale height of the thick disk 
(which is likely to lead to its contraction) would likely lower the 
fd/zo.awr ratio, irrespective of whether this additional gas is ac- 
creted by merging of gas rich satellites or through accreting gas 
from the halo or intergalactic medium. 

Also shown in Fig. [12] are the ratios for two thick disks 
formed through internal clumps scattering, which are radially 
more extended but less thick. Their ratios are on average about 2 
times larger than those of the observed thick disks. In this case, 
additional gas accretion from the galaxy halo or surroundings 
would only increase the discrepancy between observations and 
the models. 

The loss of angular momentum in the gas during the merger 
causes both the deepening of the gravitational potential and en- 
hanced star-formation in the densest regions of the disk, which 
is expected to lead to a larger scale length of the thick disk than 
of the thin disk newly for med during the mer ger This is also 
supported by observations (iPohlen et al.ll200'7l) . In order to es- 
timate the scale length of the thin disk, we consider only stars 
at disk height |z| < 1 kpc, whereas stars at greater heights 
(1 < |z| < 5 kpc) are considered to constitute the thick disk - 
as w as done in the analy sis of observations of external galaxies 
(e.g.. lNeeser et al.ll2002b . A comparison between our results and 
observations (Fig. [T3T l shows a very good overall agreement in 
the ratio of thick to thin disk scale lengths. This is a natural out- 
come of minor mergers with a primary disk containing gas and 
does not require any fine tuning. 

IStewartet~an(l2008h . in a study of the merger history of dark 
matter halos, found that virtually every MW-sized halo has ex- 
perienced at least one minor merger with m ass ratio <10:1 in 
the last 10 Gyr (see also lGuo & Whitell2008h . One minor merger 
of mass ratio <10:1, as we have shown, is not going to destroy 
the disk and leads to something that looks like a thin -n thick disk. 
However, the effect of dissipationless minor mergers on thick- 
ening the stellar disk is cumulative. If the merger frequency is 
higher, then the merger-induced disk thickening will efficiently 
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destroy the original thin disk, leading to a spheroid dominated 
early type galaxy. This predicament implies that the evolution 
history of most disk galaxies must be relatively quiet and sug- 
gests that there is a limit to how much mass can be accreted 
solely through minor mergers. However, we have yet to inves- 
tigate dissipative multiple mergers and expect that the addition 
of gas will allow some additional dissipation of the energy and 
angular momentum of the merger This would also have the ef- 
fect of erasing the perhaps unique signatures of the timing of the 
previous event an d result in a more complex evolution (see also 
iBrook et al.ll2005h . 

The question arises if metallicity or metallicity gradients 
can be used to distinguish between variou s thick disk forma- 
tion models. Within a cosmological context. ferook et al] (l2005h 
discussed that merger-induced star-formation and on-going gas 
accretion could result in a thick disk stellar population with a 
homogeneo us age and metallicity distribution. In local edge - 
on galaxies ( lMouldll2005l:ISeth et alJl2005l:lReikuba etai]|2009h . 
no-significant vertical color gradients are observed, which sug- 
gests the absence of metallicity or age gradients. Although some 
disk galaxies i n the local univ erse show strong radial metallicity 
gradients (e.g.. IShieldsll 1 9901) . they do not constitute a stringent 
test of the merger hypothesis, given the dependence on the pre- 
existing disk and the strong mixing expected. 

Our models would not necessarily predict a strong radial or 
vertical gradient in metallicity or age, for two main reasons: the 
first is that a vast majority of the stars in the thick disk at heights 
<4-5 kpc are dominated by stars originally in the thin disk, 
which implies that any radial or vertical gradients in either age or 
metallicity depend on the properties of the pre-existing disk and 
not on the properties of the merging satellite galaxy. The sec- 
ond is that mergers are effective in radial mixing, as the induced 
changes in the orbital families of the disk and in the velocity dis- 
persion would dim inish any intrins i c gradients that the origina l 
disk may have had (lHavwoodl200i:ISchonrich & Binne\^l2009h . 



Furthermore, accreting gas from either the surrounding dark 
matter halo or cosmic filaments will lower the disk scale height 
over time. This would have the impact of increasing the ratio 
of the scale lengths of the thin and thick disks. We have found 
good agreement for this ratio between our simulations and ob- 
servations of nearby galaxies (Fig. \T7\ . The agreement would 
still be reasonable if the scale height of the thick disk were to 
decrease by almost a factor of 2 (through an accretion-induced 
increase in the mass surface density of a factor of 2) in our minor 
merger simulations. For the unstable disk models, which already 
produce a ratio that is too large, a decrease in the scale height 
would make the discrepancy even worse. 

However, we have yet to investigate multiple mergers which 
are separated by episodes of significant gas accretion from the 
halo or filaments. We would expect that such gas will allow for 
additional dissipation of energy and angular momentum of any 
subsequent mergers and make the disk more robust against de- 
struction. Thus if gas accretion is important for galaxy growth, 
the limit on the number of minor mergers before the effective de- 
struction of the disk may be larger than indicated by the results 
presented here. 

The caveat is of course that the variance in the metallicity of 
the thick disk stars will have increased if the stellar population 
of the disk has evolved significantly in between merger events. 
Of course, if gas is accreted and depending on its metallicity 
and how much star-formation occurs, the metallicity of the disk 
may have increased or perhaps even decreased. Therefore, for 
an inherently stochastic process like the minor merger-induced 
thick disk formation scenario, the key factor for determining the 
characteristics of the mergers is not the me an value of the stella r 
metallicity, but its variance in abundances dRi chard et al.ll2010h . 

As observations indicate that thick disks are clearly dis- 
tinct from the thin disk co mponent (e. g..'Prochaska et al. "2OOO; 
Feltzina et al. 2003; Bensb v et ai] |2004a b; Bensby & Feltzin^ 
2006; Fuhrmann 2008) and that they must have formed early and 
relatively fast, this implies that the only viable formation scenar- 
ios are those that form thick disks early, with the thin and the 
thick disk aligned - with either a parallel or a perpendicular an- 
gular momentum vector - and with a generally rela tively larger 
scale length for the thick disks (iPohlen et al .1120071) . The minor 
merger-induced thick disk formation scenario merely requires 
a pre-existing disk in order to work. In all other scenarios the 
thick disk either relies critically on the pre-existing disk itself, 
with the thick and the thin disk at least partially co-evolving - as 
in the secular evolution models, which must also explain lagging 
and counter-rotating thick disks - or the thick disk predates the 
thin disk formation. The difference in scale lengths of the thick 
and thin disks and the weak chemical and dynamical relationship 
between the two suggest that any viable mechanism should pre- 
serve their rough alignment and have much of the mass of the 
thin disk form later Clearly, an early phase of multiple minor 
mergers would be consistent with the data if such a phase were 
to last sufficiently long to allow for some enrichment due to type 
la supemovae. Issues of timing and total accreted mass are also 
important, so as not to destroy the disk completely, which would 
then provide a natural way for the angular momentum vectors to 
remain parallel over a Hubble time. 

Because mergers can cause strong gas inflows (e.g., 
iMontuori et al.l I2OIOI) . we would also expect that the proper- 
ties of the merger-induced thick disk relate somehow to the 
bulge growth. Although not discussed in detail in this study, 
we indeed find that the bulge component does grow during 
minor mergers, thus slowly shifting the galaxies towards ear- 
lier Hubble types, and that the effects of minor mergers are 
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Fig. 13. Comparison between observations and our model re- 
sults of thick-to-thin disk scale length ratios, rd.ihick/i'd.thw, as 
function of circular velocity, iv,,., from our minor merger mod- 
els with two gas fractions, fgas = 0.1 and 0.2, and two 
clumpy gas-rich disk models (gSb-nul and gSb+u2, see § |2] 
for the notation s). Observations are shown for samples of 
late-type spirals ( Yoach im & Dalcanton" 20p"6l soUd black and 
open circles), (IWu et al.i j2002i: lAbe et aini999t iNeeser et all 
|2002|, diamond and triangles), a sample of early type galaxies 
(^Pohlen et al."2004', diamond with a cross) and the Milky Way 
(iLarsen & Hu mphrevs 2003j, square). In additio n, a cosmologi- 
cal simulation is shown (from lBrook et al.ll2004 asterisk). 



cumulative in thickening the disk. A natural outcome of mi- 
nor mergers is therefore that both the significance and scale 
height of the thick disk should be related to whether a galaxy 
is a late- or early-type disk type. Moreover, since mergers can 
have a wide range of initial orbital parameters, we would ex- 
pect to find counter-rotating thick disk stars if the interaction 
are violent enough or that the some satellites haves sufficiently 
high density to begin to contribute sig nificantly to the kinemat- 
ics of the thick disk (lOu et al.i 1201 ll) . Observations have sug- 
gested that the thick disk scale height is related to Hubble type 
dde Griis & Peletie"illl997l) . that some disk dominated late-type 
galaxies do not have significant thick disks, and that some thick 
disks counter-rotate (iMorrison et al. 1 11994 iGilmore et al.ll2002t 
lYoachim & Dalcantonll2005ll 

Finally, our galaxy models show properties typical of local 
disk galaxies. Although the details are controversial, if the MW 
formed both its ha lo and thick disk at a r ound the same time, 
over 10 Gyr ago (iProchaska et al.i l2000t iFeltzing et all 120031: 
iBensbv et al1l2004allbl : lFuhrmannll2004r then the initial condi- 
tions of the MW disk would have been probably different from 
what we assumed. However, due to the lack of knowledge of the 
properties of disks in the early universe, the initial conditions 
adopted for our galaxy models represent a reasonable starting 
point for an exploration of the impact of gas fraction, orbital pa- 
rameters and merging history on thick disk formation. 



5. Conclusions 

We used N-body numerical simulations to investigate if minor 
mergers are a plausible mechanism to form thick disks that have 
morphological and kinematic properties consistent with the ob- 
servations, in particular the vertical surface density profile and 
the disk scale height. We explored in particular the stellar excess 
formed through minor mergers. The possibility that minor merg- 
ers can scatter stars originally in the thin disk (or in the bulge) 



to large height s from the ga l axy rn id-pl ane has recent l y bee n 
pointed out by iPurcell et aU (1201 Ol) and IZolotov et aP J2010l) . 
but we have for the first time systematically studied the prop- 
erties of this excess: its vertical distribution, its dependence on 
gas fraction, orbital parameters and consecutive mergers. Our 
minor merger simulations cover a variety of initial configuration 
parameters, e.g. gas-to-stellar mass fractions {fgas = 0, 0.1, and 
0.2), initial orbital energies and orientations, two different mass 
ratio (20:1 and 10:1), as well as the cumulative effect of two 
consecutive dissipationless minor mergers. 

The main results of our study can be summarized as follows: 



- Minor interactions efficiently and gradually start heating pre- 
existing thin stellar disks from the first pericenter passage of 
the satellite galaxy, well before merging is complete. The re- 
sulting vertical surface density profile of the stellar thick disk 
is characterized by a sech function. Furthermore, a stellar ex- 
cess is visible in the profile at greater heights z > 4 - 5 kpc 
in the inner regions (r < 3Rd), which can also be fitted by an 
additional sech function with a relatively larger scale height. 
Due to this excess, a double-sech profile is necessary to fit 
the entire vertical surface density profile; 

- Usually the vertical stellar profile is dominated at all disk 
heights by stars that were initially in the disk of the primary 
galaxy. Only in a few dissipationless merger models we find 
that the stellar excess at greater heights is dominated by the 
contribution of satellite stars; 

- The excess has morphological and kinematical properties 
which are distinct generally, but see the 2x 10: 1 case in Fig. [8] 
from those of thick disk stars (however, see our most massive 
accretion simulation in Fig. [8] the 2x10:1 case); 

- The scale height of the stellar excess is constant in radius, 
while the thick disk scale height increases with radius - this 
disk flaring is a characteristic of t hick disks formed by m inor 
mergers, as shown previously bv lBoumaud et al.l(l2009l) ; 

- The scale height of the thick disk decreases with an increase 
in the gas-to-stellar mass fraction of the primary disk, while 
that of the stellar excess does not: no significant dependence 
on the gas fraction in the primary disk is found, at least for 
fgas < 0.2; 

- The effect of two consecutive mergers is cumulative, both 
on thick disk and stellar excess scale height. However, obvi- 
ously, after some number of consecutive mergers, their im- 
pact will saturate, but the exact number necessary for this 
will depend on a variety of factors such as mass ratio, orbital 
parameters, etc, which has yet to be explored; 

- Stars in the stellar excess rotate slower than stars in the thick 
disk, and their kinematics are com patible with those of high- 
a abundant stars found recently bv lNissen & Schuster! (l2010l) 
in the solar neighborhood; 

- In minor mergers the scale length of the thick disk is larger 
than that of the thin disk, consistent with observations of the 
Milky Way and other nearby galaxies. In simulations with 
higher (or lower) gas fractions the ratio of thick to thin disk 
scale lengths is lower (or higher). Whereas ratios larger than 
one are a natural outcome of minor mergers, the evolution 
of a "clumpy disk" (simulating, at least to some extent, the 
influence of secular processes on the evolution of the disk) 
produces small ratios that are generally inconsistent with ob- 
servations. Subsequent gas accretion from the halo, or along 
filaments in the inter-galactic medium, will likely decrease 
the thick disk scale height. While this does not represent a 
major problem for the minor merger scenario, for the unsta- 
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ble disk model it makes its discrepancy with observations 
worse. 

We also compared the properties of thick disks formed by 
minor mergers and through strong scattering by self-gravitating 
clumps in unstable gas rich disks (a secular "clumpy disk" sce- 
nario). We find that the clumpy disks have smaller disk scale 
heights, zo ~ 0.85 kpc, whi ch are independent o f radius (in 
agreement with the results of iBoumaud et al.ll2009l) and unlike 
the merger-driven thick disks. The thick disks formed through 
scattering do not show any stellar excess at large disk scale 
heights, unlike the merger simulations. Thus a possible way to 
distinguish between this secular thick disk formation scenario 
and minor merger models is to investigate the radial distribu- 
tion of scale heights, as well as the presence and properties of 
the excess component in the vertical stellar profile at large disk 
heights. 

Our results also put limits on the amount of accretion that 
thin disks can undergo in their lifetimes. This is not a sim- 
ple function of accreted mass, as merger-induced star formation 
within the pre-existing disk will act to both stabilize the disk 
against destruction and to limit the mass fraction of thick disk 
stars by dissipating part of the energy of the interaction. Another 
component which may play an important role in thick disk for- 
mation is the bulge component, as observations have implied that 
thick disk properties relate to the Hubble type. The existence of 
a pre-merger bulge component, or its growth during the merger 
process, could prevent the stellar disk from being tidally dis- 
torted, or weaken the disk heating, and thus result in a relatively 
less thick stellar disk in the post-merger phase. Detailed study of 
the interplay between galactic bulge growth and thick disk for- 
mation could further our knowledge of the thick disk component 
as well as galaxy evolution in general. 
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